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I . DESCRIPTION 


GRAYS is a Fortran program which will generate gray shading for the 
SD-4060 graphics device. The program produces 10 shades of gray ranging from 
. no shading at all to complete coverage of the film frame (see fig.l). The 

graphing capabilities are summarized and illustrated in figures 1 to 30. The 
• » 
figures displayed are representative of the microfilm output, but the distinc- 
tion between various intensities is much clearer on the film, expecially at the 
more intense shading. ■ 

The general idea is that given an X-coordinate, X, a width, WIDTH, a Y~ 
minumum, YMIN, and Y-maximum, YMAX , the area indicated below is shaded by one 
of the gray shades from 0 to 9. 



I— WIDTH 


(X+WIDTH,YMAX) 



YMAX-YM 


IN 


(X+WIDTH.YMIN) 


DIAGRAM 1 

II. SHADING METHOD 

The shading is accomplished by drawing one or more verticle lines 
from YMIN to YMAX in the interval determined by WIDTH. The number of lines, 
distance between lines, width of lines and concurrency of lines, are all varied 
to get the d1ffe**£ftt shades . For some shades, short horizontal lines are add- 
ed at equally spaced intervals from YMIN to YMAX as indicated above. The actual 
shading is based on a full working area of 4096 addressable points in the hori- 
zontal direction and 3072 in the vertical with a raster unit defined as the 


I 


distance between any 2 addressable points in either the vertical or horizon- 
tal direction. 

The test shading was done on a width of 15 raster units. The shading 
was then modified to give quality plots for widths varying from 11 to 21 
rasters and acceptable shading from 8 to 10 rasters. The minimum allowable 
shading width is 8 rasters, To get shading for width greater than 21 rasters 
the shading has been generalized by taking the given raster width and plotting 
adjacent blocks each of equal raster width less than or equal to 21 until the 
given width is covered. The width of the block Is the largest width 521 which 
will cover the given width in an integral number of blocks. Hence, to plot a 
given width of 22 rasters, two blocks of 11.0 rasters are plotted adjacent to 
each other while for a given width of 50 rasters, three blocks of 16.666 are 
plotted. Figures 29 and 30 illustrate this method by covering a width of 26 
rasters with 2 blocks of 13 rasters. 

III. GENERAL PROGRAM REQUIREMENTS 

The Initialization of the SD-4060 plot package is necessary as described 
in the Programmer's Reference Manual for the Integrated Graphics Software System^. 
The minimum initialization can be done by the following two statements: 

Statement!: DIMENSION AM0DE( 200) 

Statement 2: CALL M0DESG (AM0DE, ITAPE) 

Statement 1 reserves space for the 4060 MODE array which contains flags 
for various options to the 4060 plot package while Statement 2 actually does 
the initialization of the MODE array and opens the output data set, ITAPE. 

The last call after all plotting is completed is to close the output 
data set and terminate the graphics routines. 

Statement 3: CALL EXITG (AM0DE) 
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It is assumed the user is familiar with the plot routines of the 
SD-4060 Software System but it is advisable to the new or occasional 
user to review the manual for the description of the various calls. 

IV. GENERAL PROGRAM CONSTRAINTS 

A. Definitions 

Three optional statements in the SD-4060 system often included in 
programming for the 4060 plot package are listed below to define the 
variables discussed in part B of this section. 

Statement 4: CALL SETSMG {AM0DE.19.A) 

Default A = 4095. 

Statement 5: CALL OBJCTG (AM0DE,C,Y1 ,D,Y2) 

Default C = 0. 

Y1 = 0. 

D = 4095. 

Y2 = 3071 . 

Statement 6: CALL SUBJEG (AM0DE,E,Y3,F,Y4) 

Default E “ 0. 

Y3 = 0. 

F = 4095. 

Y4 = 3071 . 

Statement 4 sets the normalization factor for the plot to 4095. /A 
in the horizontal direction. Statement 5 defines the actual area of 
the scope to be used for this plot or portion of plot. Statement 6 
redefines the limits of the area of the plot determined by Statement 5 
to be in the user's coordinate system. 

B. Constraints 

The variable, WIDTH, is the width of the space to be shaded in 
subject space coordinates and should fulfill the following conditions: 
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1. WIDTH * 8. *A/4095.*(F-E)/(D-C) where A,C,D,E,F are defined 
in statements 3,4 and 5 above. 

2. The X-axis scale must be linear because of the method of shading. 

3. the intensity values must be specified to range from 0 to 9, 
with 9 corresponding to the most intense shading. If any 
value of intensity is greater than 9, the most intense shading 
is used.. 

4. WIDTH = Ax 

where Ax is the difference between 2 adjacent X co-ordinates 
and should be constant for the arrays being plotted in any 
call to GRAYS. The program allows WIDTH to vary from call to 
call, but for each call to GRAYS the width used is the same 
for all points to be plotted. If Ax is constant for all points, 
and if WIDTH is set equal to Ax, and if there are no missing 
points, then there will be complete coverage in the x-directiorr. 
If WIDTH <Ax, there will be gaps between points and if WIDTH >Ax, 
there will be overlap shading. Hence care should be taken to 
specify WIDTH = Ax. 

Shading which requires the use of varying widths poses the problem 
of having to compare intensities essentially using a different scale. 

That is, for different widths, different scales may be chosen from any 
of those in figures 1 to 28. For example, intensity 3 in figure 1 
(21 raster width) is different for intensity 3 in figure 27 (8 raster 
width). The practice of using varying raster widths to shade on the 
same graph should be used with care. It is not recommended to change 
the width for each point. 
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To find the raster width , WIDE, of the area to be shaded given 
the width in s ubject space coordinates, WIDTH* use the following 
equation: 

WIDE = WIDTH*4095./A*(D-C)/(F-E) 

Then figures 1 to 28, which illustrate the shading for various raster 

/ 

widths, can be used to get an idea of the shading which will be 
selected. 

As stated before, optimum shading is obtained for a raster width 
greater than 10, but acceptuhle shading can be obtained for raster 
widths of 8, 9 or 10. 

V. MODES OF OPERATION 

Depending on how the input arrays are dimensioned, GRAYS can 
operate in any of 3 modes. 

A. Mode 1 'i Non-Dimensioned Variables 

In this case GRAYS is called each time a point is to be plotted. 

There are no dimension statements necessary for GRAYS and the calling 
sequence is: 

CALL GRAYS ( YMIN, YMAX, WIDTH, X,INTS,1 ,AM0DE,1 ,1) 

All variables are in subject space units. 

-YMIN and YMAX are the lower and upper bounds of the interval to be shaded. 
-WIDTH is the width of the area to be shaded. 

-X is the X-coordinate of the left side of the area to be shaded. 

-INTS is the intensity value which determines the shade (0^INTS<9). 

-AM0DE is the SD-4060 mode array. 

Care should be taken to have WIDTH the same in each call to GRAYS unless 
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the user is certain it should vary (see Section IV). 
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t 
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DIAGRAM 2: 

One Call to GRAYS Plots 1 Point 


B. Mode 2: Singly-Dimensioned Arrays 

GRAYS will plot a series of points given the necessary coordinates 
YMIN (I ) » YMAX ( I ) , X(I), INTS(I) for each point' I. The calling sequence is 
DIMENSION YMIN(MR) ,YMAX(MR),X(MR),INTS(MR) 

CALL GRAYS (YMIN.YMAX, WIDTH, X, INTS, 1 ,AM0DE,M,1 ) 


All variables are in subject space coordinates. 

-YMIN ( I ) and YMAX(I) are the lower and upper bounds of the area to 
be shaded for the Ith point. 

-WIDTH is the width of the area to be shaded. 

-X(I) is the X-coordlnate of the left side of the area to be shaded for 
the Ith point, 

-INTS(I) is the intensity value which determines the shading of the Ith. 
point. ( 05 INTS (1)59) 

-AMODE is the SD-4060 mode array. 

-M is the number of points to shade in this call. (OSMSMR) 

In these definitions 1=1,..., M 

Thus it is possible to plot points anywhere in the subject space as is 
indicated by diagram 3. 
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Note that as a special case. If all X ( I ) are the same and YMIN and YMAX 
vary, a column will be plotted as In diagram 4. 
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DIAGRAM 3: DIAGRAM 4: 

One Call to GRAYS Plots M Points One Call to GRAYS to Plot a Column 


In order to plot a particular row or band, two methods could be used. 

Since to plot a band, all YMIN(I) are equal to the lower bound and all 
YMAX ( I ) are equal to the upper bound of the band, the same calling 
sequence described for this mode can be used where YM I N ( I )*YMIN( J) and 

YMAX( I )*YMAX( J) for I and J « 1 M. But since all the YMIN(I) are equal 

and all the YMAX(I) are equal, It Is not necess. ry to dimension YMAX 
or YMIN and then the calling sequence is 

DIMENSION X(MR) , INTS(MR) 

CALL GRAYS ( YMIN, YMAX .WIDTH, X , INTS .M.AMODE ,1 .MR) 

The variables are as described previously for MODE 2. Diagram 5 Illustrates 
the plotting of a single band. 


YMAX 

YMIN 



linn mi 

T I I I J 
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Diagram 5: One Call to GRAYS to Plot a Row 
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C. Mode 3: Doubly-DImeri* loned Arrays 

For efficiency In storlnq and plotting data, GRAYS Is capable of 


plotting many bands across the frame In one call. The plottlnq proceeds 
from band 1 to band 2, etc. until all bands hove been plotted. Within 
a band, the shading Is done at *he appropriate x-coordlnate. The 
plotting of a band Is Illustrated below. 



X (I.J) 

DIAGRAM 6: Plotting of Points In Band I 


The calling sequence Is 

DIMENSION YMIN(MR) .YMAX(MR) ,X(MR,MC) ,INTS(MR,MC) .NPTS(MR) 
CALL GRAYS (YMIN.YMAX, WIDTH, X , INTS .NPTS.AMODE ,MR,MC) 

-MR Is the number of bands to be plotted. 


-MC Is the maximum number of points In each band. 

-YMIN(I) Is the lower bound of band I where 1*1, MR. 

-YMAX(I) Is the upper bound of band I where 1*1, MR. 

-WIDTH Is the width of the area to be shaded (In the direction of the abscissa). 


-X ( I , J ) is the x-coordinate of the left side of the area to be shaded 
for the Jth point In band I. (The x-coordlnate of the rlqht side Is 
X(I,J)+WIDTH). 

- INTS( I , J ) Is the Intensity of the Jth point in band I. (The ranoo of 
permitted values Is 0-9). 
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-NPTS(I) Is the actual number of points to be plotted in band I 
(0<NPT$n)<MC for all I). 

VI. PROGRAM-RUN STATISTICS 

A. Timing 

IBM 360-91: The average time to plot one point Is 
CPU: 0.00094 seconds 
10: 0.00055 seconds 

Hence In 1 minute of CPU and .6 minutes of 10 approximately 63600 points 
can be plotted. 

IBM 360-75: The average time to plot one point is 
CPU: 0.003 seconds 
10: 0.0006 seconds 

On this machine in 1 minute of CPU and .2 minutes of 10 approximately 
20,000 points can be plotted. 

This is based upon plotting the 49000 points shown in figures 1 to 28. 
There is equal distribution of Intensities and widths so the actual 
timing for a particular set of data may vary slightly according to the 
raster width used and average intensity to be plotted. 

B. Size 

Complied by the Fortran H compiler with 0PT=2, GRAYS uses 4834 decimal 
bytes of core or approximately 4.8K. 

VII. PUBLICATION PRINTS 

Although 16 millimeter film is quite good for viewing plots on a micro- 
film reader, it does not produce acceptable results when an enlarged 
reproduction is made. Perhaps with sophisticated equipment which would 

lighten the dark areas and/or darken the light areas, good quality 

* » 

enlargements can be made. This has not been tried. 
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It is recommended then that 35 millimeter film be obtained of any plot 
to be reproduced. Figures 1-30 were output by the SD-4060 onto 35 nm film 
and then enlarged to 8"xl0" fine grain positive contacts using a 50 rmi 
lens at 2.6 seconds and an exposure of f/11. These were then contacted 
onto matte paper at an exposure of 3/10 sec. This method and settings 
may be used as a guide in obtaining enlarged prints. 

VIII. EXAMPLES 

A. Doubly Dimensional Arrays 

Odd numbered plots in figures 1-29: These plots were all produced 

using doubly-dimensioned arrays with 10 bands each containing 60 points 
plotted across the page. The Fortran code is outlined below for Figure 1: 
DIMENSION YMI N ( 30 ) , YMAX ( 30 ) ,X(30,100) , I NTS (30, 100) ,NPTS(30) 

DIMENSION Z(200) (amode array) 

WIDTH = 21. 

CALL M0DESG(Z,0) 

CALL SETSMG(Z,19,4095. ) 

CALL SUBJEG(Z,0.,0., 4095., 3071.) 

CALL 0BJCTG(Z,0.,0., 4095., 3071.) 

C 

C. .SET YMIN(I) AND YMAX(I). YMAX ( I ) -YMIN( I ) = 99 RASTERS 
C 

DO 1001 = 1,10 
YMIN(I) = 100. *1-50. 

100 YMAX ( I ) = YMI N ( I )+99. 

C 

C. .FILL X.INTS AND NPTS FOR 10 BANDS EACH 
C.. CONTAINING 60 POINTS 
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DO 200 J = 1,10 
NPTS(J) = 60 
NPTS(J+10) = 0 
NPTS ( J+20 ) - 0 
DO 200 I a 1 ,60 
X(J,I) = (I-1)*WIDTH 
INTS(J,I) = J-l 
200 CONTINUE 

j 

r* 

u 

CALL GRAYS(YMIN,YMAX, WIDTH, X,INTS, NPTS, Z, 30, 100) 

CALL EXITG(Z) 

STOP 

* 

END 

Note here that the dimensions were set up to handle 30 bands each containing 
a maximum of 100 points. Since only 10 bands were to be plotted NPTS< J) = 0, 

J>10. The dimensions of 30 and 100 were necessary in the call to JAYS 
in order to access the correct values. 

B. Singly Dimensional Arrays Used to Plot 3 Columns 

Figure 31 is based on a shading of 18 rasters with a scale drawn 
at the right of the grid. The scale actually is 3 points wide (3 blocks 
of 18 rasters) covering 52 rasters in the horizontal direction and each 
shade is 99 rasters high. The following code was used: j 

DIMENSION Z { 200) 

DIMENS ON VAL(20) , X ( 1 0 ) » I VAL (10) 

CALL MCDESG(Z,0) | 

CALL 0BJCTG(Z, 0.0, 0.0, 4095., 3071.) \ 

CALL SUBJEG(Z, 0.0, 0.0, 4095., 3071.) | 

[i 

, i 
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0FF=0.Q 

C.,Do 920 three times, once for each column. 

DO 920 J = 1 ,3 
DO 910 I = 1,10 
IVAL(I) = 1-1 
X(I) = 4035. +0FF 
VAL(I) = 1030.-h(I-l)*100. 

VAL(I+10) = VAL(I)+99. 

910 IF( J * EQ. 3) CALL NUMB RG(Z, 4020. ,VAL(I)+40. ,1 ,IVAL(I) ) 

CALL GRAYS (VAL(1 ) , VAL (11) ,18.0,X,IVAL,1 ,Z,10,1 ) 

920 0FF=0FF+18. 

C. Singly Dimensional Arrays to Plot Bands of Data 

The data points plotted in figure 31 were plotted with the 
following code: 

DIMENSION Z(200) 

DIMENSION EMIN ( 28 ) ,EMAX(28) ,XR(20) ,IV{20) ,YMN(4) 

DATA YMN/130. 0,855., 1580., 2305./ 

CALL M0DESG(Z,0) 

C.. REASSIGN OBJECT SPACE FOR EACH OF THE 4 PARTS OF THE PLOT 
DO 300 I =1,4 

CALL OBJCTG(Z, 200. ,YMN(I), 3900., YMN{I)+665.) 

YMIN = AL0G1 0( .4) 

CALL SUB JEG(Z, 0.0, YMIN, 210. ,3. } 

DO 250 IN = 1,28 

250 CALL G RAY S ( EM I N ( I N ) , EMAX ( IN ) ,1 .0,XR,IV ,Z,1 ,20) 

300 CONTINUE 
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IX. FIGURES 


Figures 1-30 are for display purposes of shading. The widths indicated 
are raster units per point in the horizontal direction. The microfilm 
resolution between shades is much better than is illustrated, especially 
at the higher intensities. It is seen here that the height of each point 
makes a difference at the higher intensities on these enlargements. For 
instance, a difference between the two darkest shades is distinguishable 
on the even numbered plots where the height is equal to the width, more 
so than in the odd numbered plots where the height was fixed at 99 rasters. 
On the microfilm all intensities are easily distinguishable. 

Figure 31 is a plot of actual space flight data obtained from the 
SSS-A (Explorer 45) satellite, using a raster width of 18. 


- 13 - 


X. PROGRAM LISTING OF SUBROUTINE GRAYS 

10 SUBROUTINE GRAYS(YMIM,YMAX,WIDTH,XCORD, I NTS, NPTS, Z,M, N) 

20 MMENSIOM YM I M( 1 ) , YMAX ( 1 ) , XCORD ( M, M ) , I NTS ( M, M) , Z ( 1 ) , NPTS ( 1 ) , X ( 1 4 ) 

30 INTEGER*?. I.INES(0,13)/ 

40 1 0,0, 1, 0, 0, 7 , -7,-11, 33/0,0/ 1, 1, 1, 7/ -7. +11/ 33/ 

50 2 0/0/ I, 3, 5/ 9, -a, -12, 34,0,0, 1/ 4, 5, 8, 13/+12, 34/ 

60 3 0,1, 3, 4, 6 , 1,0, -8, 21, 35,0,1, 3, 4, 6, 10, -8, 21, 35, 

70 4 0,2, 4, 6, 8, 12, +11, 25, 40,0,2, 4, C, 8, 12, +11, 24, 40, 

80 5 0,2, 4, G, 8, 12, + 11, 24, 41, 0,2, 4, 7, 9, 13, +12, + 26, 43, 

90 6 1,3, 5, 8, 7, 16,-1', +20, 52,1,3, 5, 8, 7, 16,-14,-20, 54, 

100 7 1,3, 5, 8, 7, 16,-14,-21, 56/ 

110 I MTEGER*4 l.(10)/10*0/ 

120 DATA W, IN/0,0/ 

130 0 

140 0. ..HOUSEKEEPING. .... 

150 C. ..THIS SECTION IS OOHE INITIALLY AND EACH TIME THE V/IOTH IS CHANGED 
160 C 

170 IF(W.EO.WIDTH) GO TO 300 

180 W "WIDTH 

190 C 

200 C COMPUTE WIDE IN RASTER UNITS 

210 C 

220 V/IDE»<Z(8)-Z<6>)/<Z<4)-Z(2>)*WIDTH*Z(19) 

230 imiTE(6,l) WIDTH, WIPE 

240 1 FORMAT (' GRAYS: SUBJECT SPACE WIDTH-’ , F6 .2, 

250 1 ' EQUIVALENT TO RASTER WIDTH®’ ,F6.2) 

260 C 

270 C IF WIDE GT 21 RASTERS COMPUTE NUMBER OF BLOCKS ( I PEAT) OF A 

280 C SMALLER RASTER TO BE USED IN SHADING 

290 C 

300 I PEAT-1 

310 I F(W I OE. LE. 21 . ) GO TO 150 

320 140 I PEAT » I PEAT + 1 

330 R-VMDE/IPEAT 

340 I F( R.GE. 21. ) GO TO 140 

350 WIDE-R 

360 C 

370 C ....INITIALIZE VARIABLES: 

380 C FILL ARRAY L WITH THE APPROPRIATE VALUES STORED IN ARRAY 

390 C 

400 150 DO 160 1-1,14 

410 160 X< t)®l*WIOE/15. 

420 WRI TE( 6,2) I PEAT, WIDE 

430 2 FORMAT( '+' ,73X, ' PLOTTED IN ’,13,’ BLOCK(S) OF ', 

440 1 F6 . 2, ’ RASTERS EACH.') 

450 NRAST-WI DE + . 5 
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460 I F( NRAST. LT.9) NR. AST-9 

470 NRAST«NRAST~8 

480 DO 200 1*2,10 

400 tNE»LINES( I -X, MRAST) 

500 200 L( I)«IABS( IMF.) 

510 C * 

520 C. .KEEP 4060 MODE ARRAY VALUES 

530 C... THESE WILL BE CHANG EO: THE 4060 MODE ARRAY WILL BE RESs'T UPON 

540 C... RETURN 

550 300 A30*7.( 30) 


560 A14=7.(14) 

570 A94*7.(94) 

580 A95*7.( 9 5) 

590 CALL SETSMG(7., 95,1.0) 

600 CALL S FTSMG{ 7,94,1.0) 

610 LEV-1 

620 C 

630 C... START LOOPING THROUGH ARRAYS TO BE SHAPED: 

640 C... START WITH BAND 1 (LEV-1) 

650 C... FIRST POINT IN THE BAND (1*1) 

660 C 

670 400 IMAX=NPTS( LEV) 

680 IF(N.EQ.l) IMAX-1 

690 1=1 

700 500 IFIL-0 

710 XX=XCORD( LEV, I ) 

720 NUM*M I N0 ( I NTS ( LEV, I ) +1,10) 

730 B*YMI M( LEV) 

740 T=YMAX( LEV) 

750 C 

760 C CHANGE COORDINATES FROM SUBJECT SPACE UNITS TO RASTER 

770 C UNITS 

780 CALL S F.TSMG(Z, 14,0.) 

790 CALL SCAL77.(Z,XX,B, IX, IY) 

800 B* I Y 

810 CALL SCALZZ(Z,XX,T, IX, I Y> 

820 XX«IX 

830 T= I Y 

840 CALL SETSMG(Z, 14, 1. ) 

850 CALL SF.TSMG(7.,30, .5) 

860 XF I L=XX 

870 C 


880 C... SELECT APPROPRIATE SHADING FOR GIVEN RASTER WIDTH AND INTEN .. 
890 C 

900 550 I F(VJ I D E. GE . 18 . 5) GO TO( 1400 , 40, 40, 40 , 50 , CO, 80, 70, 30 , 80) , NUM 
910 I F (WIDE. GE. 17. 5) GO TO( 1400, 40,40,40, 50,60, 30,28, 80, 80) , MUM 
920 I F(W I DE.GE, 14 , 5) GO TO ( 1400, 40,40 ,40, 50, 60 ,30 ,26 , 80, 80 ), NUM 
9 30 I F (I'M DE.GE. 12. 5) GO TO( 1400, 40, 30, 40 , 35 , 60, 80 , 25, 80, 80 ), NUM 
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940 
950 
960 
9 70 
980 
990 C 
1000 0 
1010 C 
1020 C 
1030 

1040 

1050 

1060 

1070 

1080 

1090 

1100 

1110 

1120 

1130 

1140 

1150 

1160 

1170 

1180 

1190 

1200 

1210 

1220 

1230 

1240 

1250 

1260 

1270 

1280 

1290 

1300 

1310 

1320 

1330 

1340 

1350 

1300 

1370 

1380 

1390 

1400 


I FCWIDE.GE. 11.5)00 TO( 1400, 10, 30, 2 7, 40, 28, 28, 30, 26, 80), NUM 
1 F<VM PE. GE. 10 . 5>C50 TO( 1400, 10, 30, 26, 26, 28, 30, 26, 30, 80), MOM 
I F(W I DE. GE. 9 . 5 ) GO TCK1400, 10, 25, 26, 26,28, 27, 30, 30, 80) , NUM 
I F(WI DE.GE. 8 .5) GO T0(140Q, 10, 25, 25, 27, 27, 27, 70, 30, 80) ,MUM 
RO TO (1400, 10, 25, 26, 28, 28, 28, 70, 70, 80), MUM 


. ..SHAD f NR IS DONE BY THE FOLLOWING CODE 

10 CALL SEGMTG(Z,1,XX+X(1),B,XX+X(1),T) 

XX-XX+X(8)-4. 

RO TO 30 
28 I N°l N-l 
27 I N= I N-l 
26 IN-IN-1 
25 I N° I M-l 

30 INE=(1.4*HIDE-7.0>/3. 

I F(14 . 5.GE.WI DE.AND.',,'! DE.GE. 12.5) I NE« I NF-1 
NN a (T-B)/WI OE*( INE+IABS( IN)) 

I N n 0 

I F ( NN . LT . 2 ) CALL SERMTG( Z, 1, XX+X( 7) , B + 4 . , XX+X( 7) +3 . , B+4 . ) 

I F( NN . GE . 2 ) CALL MLTPLR(Z, NN-2 , XX+X( 7) , B+l . ,XX+X( 7) +3. ,B-M . , 

1 XX+X(7),T-1.,XX+X(7)+3.,T-1.) 

I F ( MUM . EO. . 2 ) GO TO 1400 
I F( NUM. ER. 4. AND. WIDE. RE. 12. 5) GO TO 35 
I F( NUM. EQ. 7. ANO.WIDE.LT. in. 5) GO TO 65 
I F ( MUM. EQ» 6 , AMD. Wl DE. LT. 10 , 5) GO TO 37 
I F { NUM . EO . 5 . AMO .WIDE. LT .10,5) GO TO 32 

I F(Nl)M. EO,. 8 .AND. 12 . 5. GT.WI OF. .AND. W I DE.GF. 11 . 5) RO TO 40 
I F ( MUM. GE . 8 ) RO TO 70 

GO TO 40 

37 CALL MLTPLG(7,0,XX+X(6),B,XX+X(6),T,XX+X(11),R,XX+X(11),T) 

32 CALL MLTPLG( Z, 0,XX+X(4),B, XX+X( 4 ) , T, XX HI I OE, B, XX+WI OE, T) 

CALL MLTPLG(Z,L< MUM) , XX+X(7),R,XX+X(7),T,XX+X(13),R,XX+X(13),T) 
GO TO 1400 

35 CALL S EGMTG( Z, 1, XX+2 , ,B,XX+2.,T) 

40 XB=V/I DE/( 2 . *L( HUM) + 4 , ) 

XE=WI OF.+XX-XB 
XB-XX+XB 

CALL MLTPLG( Z, L( NUM) , XR, R, XB,T, XE, B, XE,T) 

GO TO 1400 

50 CALL MLTPLR(Z, L( NUM) , XX+X( 1 ) , B, XX+X( 1 ) , T, XX+X( 14 ) , B, XX+X( 14 ) , T) 
GO TO 1400 

60 CALL SEGMTG( 7, 1, XX+X( 9) , R, XX+X( 9) ,T) 

I F(WIDF. .GF.,18 .5) CALL SF.TSMG(Z, 30,1.0) 

65 CALL MLTPLG( 7, L( MUM) ,XX+X(1), R, XX+X( 1 ), T, XX+WI DE, B, XX+Wt OE, T) 

GO TO 1400 
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